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Abstract — This paper presents two new 
MAP (Maximum a Posteriori) filters for 
speckle noise reduction and a Monte Carlo pro- 
cedure for the assessment of their performance. 
In order to quantitatively evaluate the results 
obtained using these new filters, with respect to 
classical ones, a Monte Carlo extension of Lee's 
protocol is proposed. This extension of the pro- 
tocol shows that its original version leads to in- 
consistencies that hamper its use as a general 
procedure for filter assessment. Some solutions 
for these inconsistencies are proposed, and a 
consistent comparison of speckle-reducing fil- 
ters is provided. 

Keywords — Filters, image processing, im- 
age quality, simulation, speckle. 

I. INTRODUCTION 

Contemporary remote sensing relies on data from dif- 
ferent regions of the electromagnetic spectrum, the op- 
tical, infrared, and microwaves ones. Synthetic Aper- 
ture Radar (SAR) sensors are becoming more rele- 
vant in every field of research and development that 
employs remotely sensed data, since they are active 
and, thus, not requiring external sources of illumina- 
tion. They can observe the environment in a wave- 
lenght that is little or not affected at all by weather 
conditions, providing complementary information to 
the conventional optical sensors. The information 
these sensors provide is relevant for every remote sens- 
ing application, including environmental studies, an- 
thropic activities, oil spill monitoring, disaster as- 
sessment, reconnaissance, surveillance and targeting, 
among others. 



As every image obtained using coherent illumina- 
tion, as is the case of laser, sonar and ultrasound-B 
imaging, SAR images suffers from speckle noise. Such 
type of noise do not follows classical Gaussian additive, 
it is multiplicative noise. Classical techniques for noise 
reduction are thus inefficient to combat speckle see, for 
instance, Allende et al. (2001); Delignon and Pieczyn- 
ski (2002); Kuttikkad and Chellappa (2000); Medeiros 
et al. (2003) ; Touzi (2002) . 

Since speckle noise hampers the ability to identify 
objects, many techniques have been proposed to al- 
leviate this issue. Techniques are applied during the 
generation phase of the images (multilook processing, 
see Lopes et al, 1990) or after the data is available to 
the users (processing with filters). A "good" technique 
must combat speckle and, at the same time, preserve 
details as well as relevant information. 

In order to assess the performance of speckle- 
reduction techniques (multilook or filter-based), Lee 
et al. (1994) proposed a protocol. It consists of a 
phantom image corrupted by speckle noise processed 
by speckle-reduction tecniques. Measures of quality 
are computed on the images obtained, and the per- 
formance of the used tecnhique is assessed from these 
measures. This protocol can be applied to both multi- 
look or filter-based speckle-reduction procedures. In 
this paper we will discuss the use of this protocol, 
termed "Lee's protocol" , on filter-based techniques. 

This papers presents situations where Lee's proto- 
col is inadequate and should be replaced by Monte 
Carlo experiments; the outline of such simulation is 
presented. This approach aims at results that are rep- 
resentative for a collection of images, while the ones 
provided by Lee's protocol regard only one image and, 



as will be shown, can be biased and uninformative. 

Among the many approaches for speckle-reduction 
using filters, one should mention those based on the 
noise statistical properties and other general purpose 
techniques (median, mean etc.). For a comprehensive 
review of speckle filters the reader is referred to the 
works by Lee et al. (1994), by Medeiros et al. (2003) 
and to the references therein. 

Among others, Geman and Geman (1984) show the 
advantage of using stochastic models in image process- 
ing. They propose a general transformation setup: 

X — ^ Z - X 

where X represents the unobserved true data, r trans- 
formations imposed by the sensor and Z the observed 
data. The sensor typically degrades the truth by 
means of non-linear and non-invertible transforma- 
tions and with the inclusion of noise. The knowledge 
of the properties of both X and r allows building tech- 
niques for obtaining X, an estimate of X. 

With the above setup and a Bayesian approach, 
many estimators can be used to compute X as, for in- 
stance, maximum likelihood, minimum squared error, 
maximum posterior mode and maximum a posteriori 
(MAP). One of the contributions of this paper is the 
proposal of two new MAP estimators for the ground 
truth of SAR imagery, and the assessment of their per- 
formance. 

Image quality assessment in general, and fil- 
ter performance evaluation in particular, are hard 
tasks (Wang et al., 2002). Many factors are involved 
as, for instance, the true scene, the degradation and 
the type of application sought for the data. Lee's pro- 
tocol is a proposal for filter assessment based on the 
extraction of measures of quality; in its original ver- 
sion the performance of a filter is assessed using a sin- 
gle image X as input. The use of statistical models al- 
lows the proposal of Monte Carlo experiments in order 
to compute quantities that, as noise reduction perfor- 
mance, depend on many factors and are hard (or im- 
possible) to obtain directly (see, for instance, Bustos 
and Frery,1992; Robert and Casclla, 2000). We show 
here that Lee's protocol requires such a Monte Carlo 
procedure in order to precisely assess the performance 
of speckle filters. 

The paper is organized as follows. Section II 
presents the multiplicative model. In Section III two 
new speckle filters based on a Bayesian approach on 
this model are proposed, after reviewing one of the 
most widely used techniques: Lee's filter (proposed by 
Lee, 1986). Section IV presents Lee's protocol, and 
Section V describes a Monte Carlo experiment used to 
assess the new filters with respect to the one proposed 
by Lee et al. (1994); in this section we see that the 
original protocol is inadequate for filter assessment. 
Sections VI and VII present the results and conclu- 
sions. 



II. THE MULTIPLICATIVE MODEL 

Only univariate signals will be discussed here; the 
reader interested in multivariate SAR statistical mo- 
delling is referred to Freitas et al. (2005). 

Goodman (1985) provided one of the first rigo- 
rous statistical frameworks, known as "Multiplicative 
Model" for dealing with speckle noise in the context 
of laser imaging. The use of such framework has led 
to the most successful techniques for SAR data pro- 
cessing and analysis. This phenomenological model 
states that the observation in every pixel is the out- 
come of a random variable Z : fl — >• M + that, in turn, 
is the product of two independent random variables: 
X : O — > M + , the ground truth or backscatter, related 
to the intrinsic dielectric properties of the target, and 
Y : il — > M + , the speckle noise, obeying a unitary mean 
Gamma law. The distribution of the return, Z = XY, 
is completely specified by the distributions X and Y 
obey. The univariate multiplicative model began as a 
single distribution, namely the Rayleigh law, was ex- 
tended by Yueh et al. (1989) to accomodate the K 
law and later improved further by Frery et al. (1997) 
to the G distribution, that generalizes all the previous 
probability distributions. 

The density function that describes the behavior of 
the speckle noise is 

L L 

fx(x) = ^x L - 1 cxp{-Lx},L>l,x> 0, (1) 

where L is the number of looks, a parameter related 
to the visual quality of the image and that can be 
controlled to a certain extent during the generation 
of the data. An effective filter will tend to increase 
the value of this parameter, that when is estimated is 
referred to as "equivalent number of looks" . 

The most successful models for the backscatter are 
particular cases of the Generalized Inverse Gaussian 
distribution (Frery et al., 1997), being the main ones 
a constant (c) and the Gamma T(a, A), Reciprocal 
of Gamma r (a, 7) and Inverse Gaussian IG(u>,a) 
laws (for this last see Miiller et at, 2000). These mo- 
dels for the backcatter yield the following distributions 
for the return Z, respectively: 

Gamma: characterized by the density function 

/(z) = T^" z ^ lcxp{_iz/c} ' (2) 

with c, z > and L > 1, denoted T(L,L/c). 
K: whose density function is 

/co = 2 ( ^° + N (a+L)/2g °-^ 2 ^)» 

T(L)T{a) 

denoted JC(a, A, L), where z,a, A > 0, L > 1, 
and K v is the modified Bessel function of third 
kind and order v; 
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7«r(L)r(-a) (j + Lz) L - a 
where — a, 7, z > 0, L > 1, denoted £/°(a, 7, L). 
GH: characterized by the density function 
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where w, ct, 2 > 0, £ > 1, denoted G H (u, a, L). 

The use of these models in SAR image understan- 
ding has led to excellent results, as can be seen in 
the works by Mejail et at, 2003; Quartulli and Datcu, 
2004. 

New Bayesian filters for retrieving information on 
X based on the observation of outcomes of Z will be 
presented in the next section. 

III. SPECKLE FILTERS 

One of the most widely used filters for speckle reduc- 
tion is recalled, namely Lee's filter; then, the general 
setup for MAP filter is provided and two new Bayesian 
filters are derived. 

A. Lee's filter 

This filter (Lee, 1986) aims at combating either mul- 
tiplicative or additive noise or a combination of both. 
It uses the observed mean and variance in a window 
to estimate the (unobserved) true backscatter by 

x = x + b(z — x) 

where b is an estimator of the ratio of the variance of 
X to the variance of Z. If no reliable model for X is 
available, its moments have to be estimated from the 
data. The mean and the variance of the backscatter 
can be estimated in every window by x = z/y = z and 

by 

_ _ -2 2 

2 _ 2 _ Z °Y 

X ~ Z (4+1)' 

where uy is the speckle standard deviation, that can 
be easily computed by means of the equivalent number 
of looks. 

B. MAP filters 

MAP (maximum a posteriori) filters are a Bayesian 
approach to the problem of estimating the properties 
of X given the observation of Z by means of maximi- 
zing the posterior distribution of X given Z 



Since fz(z) does not depend on the sought variable x, 
the MAP estimator is defined as 

x = argmax[f z \x= x (z)fx(x)] 

or, since all quantities are positive, as 



argmax[log(./W =a: (z)) + log(/x(a:))] 



(5) 



In this work two MAP filters will be derived, those 
assuming the r _1 (a!, 7) and IG(uj,a) laws for the 
backscatter. These filters will be called MAP Q° and 
Q H . Solving equation (5) in these two situations leads 
to the following estimators: 



Lz + 7 
L + l 



a 



(6) 



and 



x = 



(L + 3/2) - y/{-L - 3/2) 2 + Aw/a{Lz + ua) 
-2uj/a 



(7) 

respectively. The equivalent number of looks L is es- 
timated beforehand for the whole image using homo- 
geneous areas (Mejail et al, 2003). The parameters 
7, a, oj and a are estimated locally using the data 
available in a small window around the position being 
filtered; among the possible estimation techniques, in 
this work we used estimators based on the moments 
of order 1/2 and 1. For maximum likelihood estima- 
tion and its numerical issues, the reader is referred to 
the work by Frery et al. (2004) ; improved inference by 
resampling is treated by Cribari-Neto et al. (2002). 

The performance of these three filters is assessed in 
the following sections. 

IV. LEE'S PROTOCOL 

Lee et.al, 1994 proposed a protocol for the performance 
assessment of speckle reduction techniques. This pro- 
tocol consists of using a phantom image (see Fi- 
gure 1(a)) corrupted by speckle noise (see Figure 1(b)) 
and obtaining measures on the filtered versions; as 
an example, Figures 1(c) and 1(d) present the results 
of two candidates for evaluation: Lee's and Q° fil- 
ters. Since the geometric properties of the uncorrupted 
phantom are known (points and strips of varying width 
from 1 to 13 pixels), it is possible to quantitatively as- 
sess the behavior of the techniques. 

The equivalent number of looks, the line contrast 
and edge preserving are the used criteria to quan- 
tify the quality of speckle reduction techniques, de- 
termined by 

Equivalent number of looks: in intensity imagery 
and homogeneous areas, it can be estimated by 
NEL = (j/az) 2 , i.e., the square of the recipro- 
cal of the coefficient of variation. 



fx\z=z(x) = 



fz\x=x{z)fx{x) 
fz(z) 



Line contrast: since the phantom has a line of one 
pixel width, the preservation of this line will be 



(a) Phantom 



(b) Corrupted phantom (c) Lee filtered data 



(d) g° filtered data 



Figure 1: Lee's Protocol phantom, speckled data and filtered images 



assessed by computing three means: those in the 
coordinates of the original line (xe) and in those 
corresponding to two lines around it (xe 1 and 
xe 2 ). The contrast is then defined as 2x£—(xe 1 + 
X£ 2 )/2, and compared with the contrast in the 
phantom image. 

Edge preserving: it is measured by means of the 
edge gradient and variance. The former is com- 
puted as the absolute difference of the means of 
strips around edges, while the latter is as the 
former but using variances instead of means. 

The best filter corresponds to the smallest value 
for all the measures, with exception of the Equivalent 
number of looks(in this case it occurs the opposit). 

It is noteworthy that the original proposal (see Lee 
et at, 1994) uses a single corrupted image to make 
the assessment. This is a twofolded limitation, first 
related to the precision and secondly to the inherent 
variability of the measures, as we illustrate with an 
example. 

A Monte Carlo version of Lee's protocol is proposed 
in next section for improved precision and significance, 
and the filters are then assessed with it. 

V. IMPROVED PROTOCOL FOR FILTER 
ASSEMENT 

The result of applying Lee's protocol to an image ob- 
tained by simulation is shown in Table 1, where the 
"best" results are highlighted in boldface. As can be 
seen, the four criteria lead to conflicting conclusions: 

• the best edge gradient is achieved by the Q° fil- 
ter but, at the same time, it has the highest va- 
riance, which is detrimental; 

• the equivalent number of looks and the edge va- 
riance point at Lee's filter as the best, but it is 
the worst in line preservation; 

• if line preservation were to be used alone, the 
best filter is Q H . 



As previously seen, Lee's protocol preconizes the use 
of quantitative measures for the assessment of speckle- 
reduction techniques. Their scale and relative impor- 
tance are not of the same order, for example the diffe- 
rence in edge variance introduced by Lee's filter and 
g H is 143.00 - 141.29 = 1.71, while the equivalent 
number of looks produced by the two MAP filters dif- 
fer by 19.04— 17.24 = 1.80, so it remains unclear how, 
if possible, to combine these results in a single conclu- 
sion. Moreover, since no distributional assumption is 
made on these measures, one does not know anything 
about the significance of seemingly different results. 



Table 1: Comparison of Speckle filters in a single image 



Speckle 


Edge 


Edge 


NEL 


Line 


filter 


Gradient 


Variance 




Pres. 


Lee 


38.46 


141.29 


25.17 


79.88 


g° 


4.11 


164.15 


19.04 


11.30 




8.97 


143.00 


17.24 


0.21 



An idea of precision of each measure, and not mere 
point estimation, becomes thus paramount in order to 
make fair filter assessments. We propose the use of a 
Monte Carlo experiment as a means for obtaining this 
relevant information. 

Such simulation experience will also help removing 
conflicting results as, for instance, the following situa- 
tion observed for the same set of parameters and two 
simulated images Z\ and Z^- Using the filtered ver- 
sions of Z\, one may conclude that Lee's filter is bet- 
ter than Q°, since the equivalent number of looks they 
produce are, respectively, 28.32 and 15.09. But if Zi 
filtered images are used, the same measures are 22.54 
and 29.88, leading to the opposite conclusion. Same 
conflicting results arise when other quality measures 
are applied, so having an idea of the variability of these 
factors is indispensably necessary. 

Since all the observed data involved have a stochas- 
tic nature it is possible to devise Monte Carlo experi- 
ments. Such in silico experiment will allow the assess- 



ment of the performance of the MAP filters provided 
by equations (6) and (7), avoiding the risk of basing 
one's decision on a single observation of a process. 

It is necessary to design a Monte Carlo experiment 
with several situations in order to provide a global as- 
sessment of filter performance. This experiment con- 
sists of simulating corrupted images as matrices of in- 
dependent samples of the Q° distribution with differ- 
ent parameters. These parameters are based on previ- 
ous experience with real data, and depict a few typi- 
cal situations often encountered when analyzing SAR 
imagery These situations stem from a constant back- 
ground (the original proposal, called "Situation 0") 
to extremely heterogeneous return (Situations 1 to 6). 
Besides simulating different types of return, various 
contrasts between the dark background and the light 
foreground were also considered. One hundred repli- 
cations were performed for each of the situations de- 
picted in Table 2. 

Table 2: Simulated situations with the Q°(a, 7) distri- 
bution 

Situation ID a 7 Background mean 

1 ^2 230 23CT 

2 -2 50 50 

3 -4 690 230 

4 -4 150 50 

5 -10 2070 230 

6 -10 450 50 



Every simulated image was subjected to Lee's and 
MAP Q° and Q H filters with windows of size 7x7, and 
the comparison was perfomed by means of the criteria 
presented in Section IV. 

VI. RESULTS 

The results obtained are summarized by means of box- 
plots. Each boxplot describes the results of one filter in 
a particular situation, using one hundred replications. 

Figure 2 shows the boxplots of the four metrics cor- 
responding to three filters in seven situations. Vertical 
axes are coded by means of the filter ('L' for Lee, 'G' 
for Q° and 'H' for Q H ) and the situation (from to 6). 
All results exhibit a notorious variability with an ex- 
ception: when both background and foreground truth 
are constant, i.e., when there is no backscatter varia- 
tion. Besides this variability, little can be said in order 
choose amongst the filters in the situations considered 
with the metrics proposed. 

In the original situation (#0), for instance, though 
there is little variability the values are almost the same 
for the four metrics. In situation #1 only edge variance 
is capable of making a discrimination, and Lee is the 
best filter (see Figure 2(d)). For situation #2, the Lee 
filter is the best one with respect to all criteria with the 



exception of edge variance; for this last criterion there 
is no difference among filters. In situation #3 only 
edge variance is capable of discrimination, pointing at 
Lee as the best filter. Regarding situation #4, Lee 
filter is the best with respect to line preservation, but 
the worst regarding edge variance. In situation #5 
both MAP filters show improvement with respect to 
Lee in both edge gradient and edge variance. Finally, 
when situation #6 is considered, though there is no 
clear evidence, Lee is better than the other two in all 
criteria but one: edge variance. 

One can conclude, then, that the original Lee's pro- 
tocol is inadequate for filter comparison when realistic 
situations are considered, e.g, when either background 
or foreground or both vary. More sensitive measures 
should be used to solve this limitation. Sant'Anna 
(1995) proposes a means to mix different metrics in a 
single scalar allowing, thus, to assessing filters with a 
grade. Wang and Bovik (2002) also propose a scalar 
metric for image quality assessment. 

VII. CONCLUSIONS 

This paper presented two new MAP filters for speckle 
noise reduction and a Monte Carlo experiment that 
improves the original proposal. When trying to assess 
their perfomance using this simulation setup, we no- 
ticed that this procedure is inadequate to deal with 
situations more realistic than constant background. 

Previous works (see, for instance, Medeiros et al. 
(2003); Sant'Anna (1995) relate the superior behavior 
of MAP filters with respect to Lee's filter, but the high 
variability of the measures proposed by Lee's protocol 
hampers a straightforward comparison. 

Other measures are showing a good potential for this 
assessment, such as the image index quality provided 
by Oliver and Quegan (1998) and the edge preserving 
index proposed by Sattar et al. (1997). 

Monte Carlo experiments must be devised within 
the framework of statistical modelling in order to make 
sensible comparisons among filters. 
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